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ABSTRACT 


One of the more significant technological problems associated with 
the orbital operation of large astronomical telescopes is the fabrication 
and maintenance of the primary mirror surface to the tolerance required 
for diffraction-limited performance. An interesting approach to the 
solution of this problem involves continuously measuring and automa- 
tically correcting the optical surface of a thin defonnable mirror by 

i i 

means of discrete actuators located on its rear surface. The real- 
ization of diffraction-limited performance from a telescope in space 

t ’ ’ 

by this method rests on the ability of the designer to achieve extremely 

accurate control of a highly complex^ interacting, multivariable system. 

This paper presents the results of a detailed study of the discrete 

control of linear distributed systems with specific application to 

the design of a practical controller for a plant representative of a 

telescope primary mirror for an orbiting astronomical observatory. 

The problem of controlling the distributed plant is treated by 

employing modal techniques to represent variations in the optical 

figure. Distortion of the mirror surface, which arises primarily from 

thermal gradients, is countered by actuators working against a backing 

structure to apply a corrective force distribution to the controlled 

surface. Each displacement actuator is in series with a spring attached 

ii 



to th.e mirror by means of a pad, intentionally introduced to restrict 

tbe excitation of high-order modes. Control is then exerted over 

a finite number (equal to the number of actuators) of the most 

significant modes. ’ ' , 

ThroTigh the application of the modal expansion technique the mirror 

equation of motion is transformed to a set of uncoiipled^ linear, time- 

invariant, ordinary differential equations. Kie desired dynamic 

response and static accuracy may then be achieved by the application of 

classical single-variable design techniques. Th6 formulation of a ' 

quadratic performance index which Incorporates a meastiM of unage quality 

permits detennination of the trade-off between the- number of actuators 

1 • , 

and optical purity. A criterion for defining actuator placement and 
pad size is presented which minimizes the tendency of the controller to 


excite the unmonitored modes. 
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I. imRODUCTION 


Tlie class of distributed plants considered in tMs paper is 
restricted to those described by linear, time-invariant, separable 
operators where control is derived from a finite number of discrete 
inputs. Application of the modal expansion approach(^) converts the 
distributed-parameter problem to one of a multivariable nature which 
readily .yields to decoupling techniques.'’ * Classical single-variable 
control methods are applied to decoupled system dynamics defined in 
terms of the eigenvalues of the linear operator whose eigenfunctions 
are assumed to form a complete O3rthonormal set. While the results 
are applicable to the general problem of controlling linear distributed 
systems, the motivation for the study is a direct resiilt of its 
relevance to one of the central problems of large orbiting telescope 
technology . 

Elimination of the effects of the earth’s atmosphere give orbiting 
telescopes significant advantages relative to earth-based telescopes. 
Figure 1 displays percent transmission of the incident radiation 
through the earth's atmosphere as a function of wavelength. In 
the portion of the spectrum shown, the earth’s atmosphere is opaque 

to radiation shorter than 5 X 10"^ cm and longer than 5 X 10“^ cm, -with 

■? 

numerous gaps between these extremes. Since the entire .spectrum of 

if ' 

radiated energy becomes available for study in an atmosphere-free 

t 

environment the additional spectral coverage would permit studies 

' ' (^) 

involving galactic nebulae and cool stars ( stellar evolution) . 



TRANSMISSION % 
lOOr 



WAVELENGTH- cm 


Figure 1.- Transmission of incident radiation through the earth's atmosphere. 
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Secondly, refraction anomalies associated with the turbulence of the 

earth's atmosphere limit the resolving power of earth-based telescopes; 

consequently, the 200-inch (508 cm) Hale telescope at Mt. Palomar has 

no better resolution than a high quality telescope of approximately 

15 inches (58 cm) diameter. ' Placing a 120- inch (5 meter) telescope 

in orbit would yield an increase in resolution of at least a factor of 

7 relative to land-based telescopes and 5 relative to present space 

telescopes. This significant improvement is useful in studying double 

{^) 

stars and in planetary photography.' ' Figure 2 shows two views of 
the same portion of the lunar surface. The photograph on the right 
was teiken by the 120-inch (3 meter) reflector at the Lick Observatory 
and represents about the best resolution available from an earth-based 
telescope. On the left is a Lunar Orbiter photograph of the same area 
which has been selected to show the resolution available from the same 
size telescope outside of the atmosphere. Further, elimination of the 
background glow associated with the atmosphere would permit longer 
exposxire times in celestial photography, effectively enabling astronomers 
to see deeper into space. Realization of these advantages is contingent 
upon solving several technological problems. 

One of the more significant of these problems and the one which 
motivated the present investigation is the fabrication and maintenance 
of the primary mirror surface to the tolerance required for diffraction- 
limited performance. For the purpose of this paper, diffraction-limited 
performance is defined as being obtained when the rms figure error 
over the mirror surface is less than a fiftieth of a wavelength, which 
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Figure 2.- Two views of the l\inar surface indicating the increased 
resolution available from an orbiting telescope. 
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at 6528 A is l/2-|i inch (1.27 X 10”^ cm). It is extremely difficult 
to achieve this accuracy with the monolithic mirror normally used in 
telescope applications as uncorrected thermal gradients, fabrication in 
a one "g" and operation in a zero "g" environment, and materiail reaction 
to stresses introduced during figuring all tend to cause distortion of 
the mirror surface. 

An interesting approach which has been investigated as a means of 
solving the problem involves measuring the deviation of the mirror 
figure from the desired shape, generating the necessary control signals, 
and applying these signals to pttysically align the mirror to the desired 
shape. This concept, which is shown schematically in Fig. 5^ bas been 
investigated with both a segmented and a thin deformable mirror. The 
segmented mirror consists of a nvimber of individual pieces or segments. 
This approach was selected because many of the effects causing surface 
deformations are reduced on mirrors of small size. In a Cassegrain 
telescope the incoming light is reflected by the primary £ind secondary 
mirror to focus behind the primary. In this case, the pr ima ry is also 
illuminated by monochromatic light which is returned to the mirror 
figure error sensor. The error signal is then processed and applied 
to the actuators which correct each segment in tilt and focus. 

This concept has been successfully applied to a small segmented 
mirror. However, the construction of large segmented primary 
mirrors from individually fabricated off-axis portions of a paraboloid 
(matched in focal length) would require new techniques of fabricating 


Primary mirror 
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Control 


electronics 


^ ^ '[^ming /fgh^ 


Secondary mirror 


Mirror figure- 
error sensor 
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Figure 5*- Schematic representation of a system for controlling the 
optical surface of a telescope primary mirror. 
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diffraction-limited aspherics. In addition, the segmented mirror 
consists of a nimiber of monolithic mirrors each subject to the same 
limitations as a monolithic mirror. While these limitations are 
reduced in the relatively small segments, they are still of stiff icient 
magnitude to require that the segments be fabricated from a substance 
which exhibits a very high degree of material stability. This is 
because the segments are corrected in tilt and focus only and any 
warping of a segment cannot be completely corrected. 

An sd-temate approach involves a thin, continuous -surface, 
deformable mirror that can be stressed into the desired shape by means 
of a large number of actuators arrayed across its rear surface. The 
greater control flexibility inherent in this approach shows promise 
of relaxing the material stability and fabrication tolerance require- 
ments. A laboratory model of a thin mirror is shown in Pig. 4 alongside 
a conventional monolithic mirror. The thin mirror is 30 inches (?6.2 cm) 
in diameter and 0.5 inches (1.27 cm) thick. Actuators were located 
every 5-5 A inches (9*^ cm) over the rear surface of the mirror and 

/ ON 

were used to apply a corrective force distribution.' ' Preliminary 
operation of the system indicates that the actuators were able to 
reduce the initial figuring error to acceptable levels. This is 
illustrated in Fig. 5 which displays photographs of the diffraction 
pattern of the mirror before and during control. 

While this application appears to demonstrate the ability of the 
thin mirror approach to enhance the telescope performance, a number of 
areas exist in which an inqiroved design technique would be of significant 



Figure 4 .- Comparison of a thin defonnable mirror with a mirror of 
normal thickness-to-diameter ratio. 




Figure Comparison of the diffraction pattern of the uncontrolled 
mirror with that of the mirror d\iring control. 


VO 
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value. In present design efforts the dynamics of the plant have been 
largely ignored. This is a significant factor in the control of larger 
mirrors which have very li^tly damped low frequency resonances. 
Selection and placement of the actuators Is presently done on an 
arbxtrary basis as a result of the present limited .ability to relate 
the effects of these design decisions to. system performance, in • . 

additionj the ability to more completely incorporate information on 

^ 1 

the disturbance characteristics, to the extent that it becomes^ 
available, is desirable. 

The p\irpose of this paper is to present a general theory for the 

1 ( 

discrete control of a distributed-pa-rameter system and extend the modal 
expansion technique to completely specify system performance. 

In the past little research effort has been devoted to the problem of 
obtaining discrete control of distributed-parameter systems and the 
results which have been obtained -by Gould and Murry-Lasso^^^ are limited 
to plants which have finite modal content. In this thesis the entire 
modal structure of the plant is considered for the problem of obtaining 
a specified level of performance while minimizing the number of control 
inputs to the plant. The method is demonstrated with examsples and 
results are presented for a plant representative of a thin deformable 
mirror. 



II. MODAL EXPAHSION !EEC!IP0IQI1E 


A schematic Tepresentation of the plant imder consideration. is 

, M- 

t 

shown in Fig. 6 where L is a linear^ time -invariant, separable operator. 
A typical example of a distributed system is the thin rectangular plate* 
of Fig. 7 , where the deflection normal to the midplane of the plate 
w(x,y,t) results from the application of a transverse load density 
p(x,y,t). The equation of motion of the plate is given by^^^ 


S V^w(x,y,t) + p — w(x,y,t) = p(x,y,t) 
&b2 



where is the Laplacian in Cartesian coordinates, p the mass per unit 
area, and S = yh5/l2(l - v^) is the flexural stiffness modulus of the 
plate. The deflection, w(x,y,t), is assmed separable, i.e., 

w(x,y,t) = c^(t)u^(x,y) . H 

Substituting into the homogeneous form of H yields, for a 
homogeneous, plate of uniform thickness. 



+ co.Cj (t) = 0 



and 
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p 

where is a constant whose value is specified hy the houndary 

’ 1 

conditions. For the simply supported plate, which has houndary 
conditions^ ■ < ' 


and 


■u-i(O^y) - u^(a,y) = Uj^Cx^O) = u^(x,h) 


= O', 





u. 

1 


+ V 



X 


= 0 , 


0,a 




where' m, n = l, 2, ... are the mode indices corresponding to the 
eigenfunctions 


^i = ^ (sin m J x)(sin n | y) [pb] 

whxch form a complete orthonormal set (Appendix A) . Consequently, 
the general solution to Eq. W is 





(x,y,t) c^(t)u^(x,y) = ^ Cj^(t)Ujjj^(x,y) ; : [6^ 

m,n=I 


i=l 


and since the form a complete set the loading may he exp^ded in a 
imifonnly convergent series of the form 


W 

P(x,y,t) =^a^(t)uj^(x,y) = ^ a^(t)Ujjjj^(x,y) 


i=l 


m,n=l 


where 


>i(t) -J'f 'ly 


[6c] 


and r is the spatial region in which the plant is defined. 

Substituting and into jjL^ and taking the Laplace 

transform with respect to time yields 


^ (s^ + CD^)c^(s)u^(x;,y) = a^(s)u^.(x^Y) ‘[^73 

i=l i=l 

where s is the Laplacian operator andj for convenience, the same symbol 
is used to denote a variable and its Laplace transform. Since the 
Uj^(x,y) are independent the coefficients may be equated yielding 


Cj^(s) = 



1/p 




\(s) = 



15 


Based on 'Eqs. Q6a-c^ and the plant shown in Fig. 6 may he 

redrawn as shown in Pig. 8 which is a modal expansion of the distributed 
plant • 



Figure 8.- Modal representation ofs the plant. . 

> 

More generally^, Fig. 6 represents the functional equation ' ' - ' 

4, ^ *■ 

^ * I 

w(z,t) = Lp(z,t) ' ! > ' ’ H 

where z represents a general spatial variable (in one or more 
dimensions) and the operator L operates on functions of time and 
distance. Laplace transforming Eq.. ^9^ "with respect to time yields 

w(zjs) - Lp(zjs) [10] 

where the eigenvalues of L satisfy 


Lu^(z) = Aj^(s)UjL(z). 



Since the u^(z) are assumed to form a complete set, the separable 
functions w(z,s) and p(z,s) can be expanded as 







l6 


w(z,s) 


GO 

= ^ Ci(s)Ui(z) 


i=l 



and 


Substituting 


00 



a^(s)u^(z). 
into Qlo] yields 





i=L 

and consequently. 


i=l 


i=l 



c^(s) = \(s)a^(s) 



with a^(s) defined by the transfom of Eq.. Figure 8 is then 

the general modal representation of the -class of distributed plants 

under investigation. In a function space where the eigenfunctions of 

> 

L are used as the coordinates the system is represented by the infinite 
diagonal matrix * ” 


A(s) 





f 


In this reference frame, control of the plant output can bte readily 

, ' t , ‘ ' 

achieved by individually compensating each element of the diagonal 
matrix, as shown in Fig. 9. 
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error which can he reduced hy appropriate control one mode at a time 

as indicated in Pig. 9 • Relating the original signals in the system 

to those of the decoupled reference frame ^ the control system structure 

becomes that shown in Pig. 10 where e and a are column matrices whose 

elements are the modal coefficients e^_ and aj_, respectively. In 

practice the situation illustrated in Fig. 10 can only be approximated. 

The function of the analyzer is to determine the modal content of the 

optical surface error. The decoupled controller dynamicsj represented 

♦ 

by the matrix d(s) = diag [^d^(s)^, is determined on the basis of 

> 

standard design techniques (see Pig. 9 ) to achieve a satisfactory 

performance level. For a 'well-ground mirror the need for corrective 

action diminishes as the mode 'number^ i, increases and control can be 

> 1 

\ 

reasonably be restricted to the significant modes.* The N controlled 

modes are denoted by the output and error N vectors and e^ in Pig. 11. 

The finite (n x R) controller matrix, is represented by D*^. ' The 

function of the load synthesizer is to place an. appropriate force 

. * • * w 

distribution on the plate to correct for the modal errors in e'^. Since 
the remaining modes are unmonitored (no corrective action taken), the 
ideal force distribution applied by the loading mechanism is 

N 

i=l 

with none of the uncontrolled modes excited. Physically this corrective 
loading is applied by finite number of control manipulators which, in 




Figure lO-.r Controller representation in the original reference frame. 


w (z,s) 
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the mirror problem^ are comprised of displacement actuators in senes 
with a spring acting against a hacking structure. The spring is attached 
to the mirror by means of a pad intentionally introduced to restrict 
the excitation of uncontrolled modes (a point discussed in detail later 
in the paper). The actual (non- ideal) force density applied by these 
N actuators is given by 


IT IT 

p(z^t) Pj(z,t) aj(t)3j(z) [ 17 ] 

.3=1 J=1 

where p^ represents the force distribution resulting from the 
actuator^ and the last step results under the assun^tion that each 
applied force distribution is separable in time and' distance. Expanding 
each of the PjC^) in terms of the eigenfxmctions^ Uj_(z)^ Eq. [XT] becomes 


where 





and comparison with reveals 



[l8b] 


[ 19 a] 
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In matrix form this relation between the mode force' coefficients and 
the actuator signals becomes 


a 


Erf 


,N 



where H is an <» x N matrix and is an IT vector. Equation [I9b] 
indicates the control elements excite all modes. Since only H of the 
more significant modes are controlled, Eg.. Q-90 is partitioned as 


'a®' 


>■ 

a's 


■# 



— — 





where is an If X N matrix, an <» x If matrix, a^ an N vector 

•D 

corresponding to the controlled modes, and a'^ accomts for the remaining 

modal force coefficients. To provide the desired corrective vector 
If 

a the actuator inputs are given by 


- [ h “]-1 a ® [ 21 ] 

where it is assumed that the actuator locations insure is non- 
singular - this point is discussed in detail in a subseq.uent section. 
Partitioning the matrix representing the plant dynamics into components 

corresponding to the controlled and uncontrolled modes, the overall 

N R 

system becomes that shown in Fig. 11 where disturbances q and q acting 
on the plant are included as equivalent displacements. When the mode 
number, i, is ordered with increasing frequency of vibration, ui = 
the plant inherently performs modal filtering which attenuates the 
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higher inodes so 'that their contribution to the' errof (mirror distortion) 
rapLdly becomes negligible as the mode number increases. Consequently, 
it is often assvtmed that only the first W modes are present,' i-e. , 
c^ and are Identically zero,^ and considerable simplification results. 
For example , let be an K vector defined by the output at’ N different 
points. That is, 



N 

w 

= col w(z., 
J 

t) 

where z. represents a measurement point. 
0 

In terms of 

displacement coefficients 





w^ = uV 


where 

r 




u^(z^) 

UgC^i) . 


= ■ 

m 

• 

• 

* 


♦ 

^i(%) 

♦ 

• 

* • ^]j( 




25a 


[25b] 


Under these conditions the mode analyzer becomes simply ,an operation 
on the K measurements; specifically 

-1 


W 






w 


where the sensors are located at positions to insure is nonsinguiar. 
The control structure of Fig. 11 then reduces to the H x N multivariable 
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system^^) shown in Fig. 12. This idealized representation^ valid when 
the effects of the higher modes can be safely neglected^ was derived 
"by Gould and Murray-Lasso and is treated in detail in reference 1. . 

Control of the low-order high-amplitude inodes as indicated in Pxg. 12^ 

' • ( 13 ) 

or in the decoupled form of Pig. '9j p3resents the classic problem' 

of controlling a resonant plant ‘With a, limited control effort (restricted 

actuator throw). While for large primary mirrors with low resonant 

f V . ^ ' 

frequencies this may be a sub'stantialr problem^ 'in the present paper 

i I ♦ 

it is assumed that the disturbances q.^(t) are slowly varying and of 
sufficientljr small amplitude that any desired degree of control can 
be realized. Negligible contribution from the uncontrolled 'modes can 
generally be assured by permitting N to be arbitrarily large. However^ 
in the present problem extremely accurate control of the optical surface 
of the thin deformable mirror for diffraction-limited performance is 
required with a minimum number of actuators. Under these conditions . 
the effects of e^ and a^ are not negligible ^ but, in fact, represent 
the' most significant system errors and the most important factors in 
evaluating design tradeoffs. 




Figure 12. - Idealized N X H multivariable system resulting from 
a finite modal expansion of the distributed system. 
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III. ACTUATOR PAD SIZE AMD LOCATION AS DESIGN FACTORS 


The ^controlled modes enter the problem in two major ways. First, 
the actuators excite not only the controlled modes but, in general, 
all modes. To demonstrate, first separate the image quality index 
into two parts, i.e., 

■'l - \ * ""e [=50 


where 


M 



i=l 


accounts for the error in the controlled modes and 


GO 



i=M+l 


accounts for the remnant error of imcontrolled modes. As larger 
actuator displacements are commanded in order to reduce J^j. to smaller 
and smaller values Jp increases due to the effect of a^. Secondly, . 
unless direct measurements of the modes are made, a limitation on the 
ability of the displacement sensors to obtain, an uncorrupted estimate 
of the N controlled modes results from the presence of e . t 

Actuator Pad Size 

The function of the actuators, as indicated in the previous 

section, is to apply forces to reduce e^ while minimizing the 

excitation of the uncontrolled modes, i.e., ideally 

26 
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H = 



nonsingular. 



In the mirror problem the control manipulators are modeled by a 
displacement actuator working against a backing plate and a relatively 
soft spring which is attached to the mirror by means of a pad. The 
ability to approach the situation on Eq.. [^2^ is governed by pad size 
and location, which are factors under the influence of the designer. 

To illustrate the effect of pad size, consider the rectangular plate 
of Fig. 7 . The pads are assumed to be rectangular in shape and 
located as shown in Fig. 15* 







Figure 15 •- Pad. shape, size, and location. 
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The amounh of force appHed "by fhe actuator is determined by the 
product of the spring constant K and displacement actuator position. 

This force loads the plate as Indicated by Eq,. [17] with 

Kf^Cx^y) 

3j(x,y) = 

f^(x,y)dx dy 




where f^(xjy) is the distribution of the force and aj(t) is the control 
input. Considering a force distribution that is constant over the area 
of the pad. Eg,. [l8b] yields 




K 




: / / z 


nwx . n:ity . , 

sin sin — — dx dy 

ab a b 


pad 

area 


4k . “‘“"j . 

— sin sin 

ab a 



mrt(Ax) A 

sin — r ^ 

2a 

n^ty. 

b 

—(Ac). 
L 2a 


sin 


2b 


not 


2b 





where i indexes the mode m,n. For this special case with constant pad 
size (i.e., (Ax)^ = (Ax)j = and (Ay)^ = “lay 

be rewritten in matrix form as 



H = GU» 
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where U is the HT x » matrix, 

Ui(z^) UgCz^),, 


U = 




Ui(Zjj) UgCZjj) . . . . 


the prime denotes the transpose of a matrix^ denotes the ^ point 

* ? 


and 


G = diag ^ 



gi = Sum = K 


jnirA 

a 

sin 

2a 


• miiA 
a 

2a 


n^- 

Sin 

2b 




Figure l4 contains a plot of (sin |)/|. Assuming the controlled modes 


“ - “max - “max dimensions are 5; 


< 


and 


Ab A 




“max 


■, the maximum value of the argument for one of the controlled 


modes is | = it/2 which occurs when adjacent pads touch. The attenuation 


of the higher order modes by the pad is apparent from this diagram as 
the elements of G decrease rapidly for i > N and, in turn, decrease 
the output levels of approaching the idealized condition of Eq.. [26], 



50 


s?n(& 

i 



Figure l4. 


Plot of sin (|)/S illustrating the filtering action of pad. 


She prefilter action of the pad is complemented by the transmission 
properties of the plate itself. For the rectangular plate of Fig» 6 
the relation between the applied loading and displacement output for 
the i^^^ mode is 


Ci(s) 


a^(s)/p 

~2 2 
s + 6 s + -o:£ 



where a small amount of damping has heen included. In response to a 
step input the steady transmission is 


lim fi 
t CO a^ 
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where oi? is given by Eg.. M- Assuming a 0.5 inch (1.27 cm) thick, 
30-inch ( 76.2 cm) square plate with a Young’s modulus of 10 pounds 
per square inch (70.3 x 10 gm/cm^) and a Poisson’s ratio of 0.2, 
this factor is 


lim °i _ 7.73 X 10“2 

t a. / 2 , 2x2 

X (m + n ) 



A byproduct of the pad's desired effect on ^ is a decrease in the 
transmission properties of H^. As a consequence, an increased effort 
is required to deflect the surface. This is readily demonstrated for 
the simply supported rectangular plate where 

. sV [ 34 ] 

and with appropriate actuator placement U*^ is orthonormal- (see 
Appendix B) . Consider the expected' value of the norm of the actuator 
displacement vector given by 


a a 




-1 W 
a 




i'j t 

Thus, as the elements of G are decreased, the required control 
displacement and force increases. Since/ the plate^ Itself was Shown 
to perform substantial filtering/ a, compromise on the final pad shape 
and size is normally employed-. 
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Actuator Location 

*R * 

Another critical factor on the H and H , matrices Is the placement 
of the actuators. The actuator locations must he selected to prevent 
the occurrence of a singular matrix. Since the actuator locations 
most directly affect the matrix, Icnowledge of the mode shapes will 
generally permit the selection of the location of the actioators. This 
situation is detailed in the examples. Evaluation of the determinant 
may he used to verify the invertlhility of H*^. In addition, if a row 
of is zero, the mode corresponding to that row cannot he excited. 

This is approximately achieved for most plants hy locating pads at the 
zeros (nodes) of the mode in question. For the rectangular plate, this 
is exactly achieved since the influence of pad location of H is 
delineated hy U’, see Eq.. . When the modes are ranked in order of 
importance, the desired goal is to control the first IT and null the 
next highest modes; however, this is not usually possible and trade- 
offs are required. For example, it may he possible to control a set 
of modes that are not the N most significant hut he ahle to preclude 
excitation of the next highest modes or, alternatively, control the IT 
most significant modes hut not have the ability to preclude excitation 
of the most dominant uncontrolled modes. Since the controlled modes may 
he reduced to any desired level at the cost of some increase in the 
amplitudes of the unmonitored modes, the/ contribution of these higher 
modes to system error represents the most critical factor In system 
design. Determination of the trade-offs in actuator size and placement 



is obtained through evaluation of the system perfonaance Index, Jj, 
which is discussed helow and illustrated later in two examples. 



IV. STAJIG ESRTORMAHCE 


The design objective is to minimize the image index defined as the 
expectation of the integral square surface error. As shown in Eq. 
this is a function of the mean square values of the modal error 
■coefficients. Determination of Jj as 
on the nature, particularly the spectral content, of the disturbances. 
For the application to the control of a deformable primary mirror of 
an orbiting telescope, it is anticipated that the primary error sources 
will be initial figuring errors and relatively slowly time -varying 
thermal gradients. In this context it is reasonable to expect that tie 
system will generally be performing at or near its static values. 

In Fig. 11 the surface deformation due to the disturbances is 
defined in terms of its displacement modal expansion coefficients, q. 

No loss in generality results from considering the disturbances to be 
displacements since equivalent force distributions could be assumed. 
With reference to Fig. 11, and with oij,(z,s) zero, the error in the 
control modes is given by 

e” = - [IH- aVJ"- g”. [36] 

For the static situation this reduces to 


given by Eqs. |^25a-(0 is depeident 


-1 


1 + 


% 


0 < i < N 



34 
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where Kj^ is the loop gain (type 0 system) for the i^^ mode. When, the 
loop transmission contains a pure integration (type 1 system), -> «>. 
Thus the" error in the controlled modes can theoretically he reduced to 
any arbitrarily small value. The expected value for Jjj then becomes 





where is the rms value of the static (or slowly vaiying) disturbance 
The error in the xincontrolled modes is given by 


e« = - 4 ® - aW 

= - 4® + a%'>’[e'G'VCi + aWJV [39] 

which 'under static conditions becomes 




1 + Ki 


IT + 1 ^ i 


where ^ is an element of the matrix 

I.. Id j ' 


t 





’ AssWing'.the modal coefficients of the disturbance are uncorrelated, i.e 




E(qj-qj) = 


1 7^ d 
i = d 



36 


the expected error in the uncontrolled >modes hecomes 



where is the value of the disturbance error in the higher modes 
without control. The second term. is clearly positive and represents 
the increase in that results from the actuator displacements 
required to control the errors in the first U modes. Since is 
finite, the series converges and the order of summation may be reversed, 
yielding 



where 

CO 

j=W+l 

is a constant dependent on the design factors of actuator pad size and 
location as well as the natural modal filtering performed by the jplant. 
•Combining Eq.s. [25a-c[] , [ 58 ^ ^ and Ql-33 yields as the imaging index 



Tfcie first term is unaffected by the control action, the second decreases 
• as a result of the control action and the third increases. 
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^Since is constant ^ minimization of Eg.* M with respect to 

i Ki yields optimuni' perfomance when 
‘ 




and 


IT 


U 


“-opt 


* Z (r^ “I - X 

i=l ' ^ 1=1 


. T 1 + (pf 
1=1 ^1 


2 

2 ^ , 
% + 


[46] 


with the controlled and uncontrolled components given by 




opt 


N IT 2 

[E (mi) \ ‘ E(ii^ 

i— 1 i=l ^ 


i/ "I 


M 


and 


N IT 2 

J-R =y ^ a§ + Jr =y tj2 + Jr . ptsl 

^(1 + K^)2 “^-i ^(1 + q>2)2 1i •- -* 


i=l 
IT 


Recalling that without control ~ ^ that for large 

° i=l 


loop gain the error in each term of JR^p^ "to control is approximately 
1 /Kj^ while each term in is reduced by an additional 

ThuS; if significant improvement is to be gained in the optical surface 
by the aboye method, IT must be selected sufficiently large and the 
actuator . size and! placement such that cp? = l/k. « 1. In the mirror 
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IHie required value of N to achieve the desired rms figuring error can 
he minimized by the selection of pad size and location whose effect 
is manifest through the parameters cp|. 



V. MOI® ESTBamOW ERROR 


The preceding analysis assumed ideal measurement of the controlled 
modal variables however, in many applications it is neither 
practical nor possible to obtain direct measurements. In these cases 
an estimate is often derived from a spatial sampling of the distributed 
output. This is the case in the mirror problem where the most commonly 
used measurement of the optical surface is performed by the interference 
method illustrated in Fig. 15 . This mirror figure error sensor is 
a modified foim of a Twyman-Green interferometer. In this interferometer 
two plane wavefront beams are formed from a common coherent source. 

One beam is reflected from a reference flat while the second is con- 
verted to a s^erical wavefront whose center of curvature is that of 
the mirror. This wavefront is returned by the mirror and forms an 
interference pattern with the reference beam which is focused on the 
n discrete individual sensors. Periodic motion of the reference flat 
produces a sine wave of identical frequency at each detector. This 
converts the error determination from an amplitude to a phase measure- 
ment and permits the required sensitivity to 'be achieved. - ' 

Under conditions where no modes except the first 31J exist, the 
relation of the modal coefficients to the N measured values is given 
by Eq. which for the mode error is' 


•"meas 




.-1 


w= 


R 
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Beam expander 
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Movable reference 
flat" (phase shifter) 


Primary mirror 


Figure 15.“ Schematic of the mirror figure error sensor. 
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However, the presence of the hi^-order inodes deteriorates this 
measurement since the actual measured vector is defined by 

+ A® [51] 

where is the H x » matrix 




Ideal 








and substituting into and -rearranging yields ' , * 

a" = - [l + aVJV . [l + aVJ\W. •• [55] 

For a type 0 system tinder static conditions,. Eg.. [p5^ becomes-, 


e. 

1 


<31 Ki ^ 

1 + “ 1 + i 


i ^ K 



or, for a type 1 system 

e® = - €» = . [#]"VaS. [set] 

Attention is now given to tbe error defined by Eq.. since, as 

previously established, most practical systems would possess an 
infinite loop gain. Two distinct cases are now considered: one -where 

is dominated by the disturbances acting on the plant, i.e., 





and, secondly, 'where the error in the higher-order modes contributed 
by the disturbance vector, is negligible but the error introduced 
by the control effort is significant, i.e., 

= - A%^o?^ = -1- 

= (I + ilr|y3 [ 58 ] 

For ease in later calculations it is assumed that iP is negligible 

compared to I and that is nearly orthogonal permitting the contributica 
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to Jjj arising from to te determined from the noim of the measurement 
error at the W measurement points That is^ 


•where 


% 

meas 


= , eQI'4] 




case I 

*' » f- 

case II 


If the modal coefficients of the disturbances are' uncorrelated then 
. . ■> '» 

E<1* |39§] case I become B ’ ^ , 


rn 




i=N+l! 


n=l 





M 

The sensor locations are chosen to minimize 

the more significant hi^ order modes as possible. This result is ' 
consistent -with that concluded with regard to actuator placement and^ 
consequently, the criteria for actuator and sensor placement are 
identical. Since increasing the number of sensors is relatively 
inexpensive, in many cases it will he desirable to have more sensors 
than actuators (or controlled modes) . If M > N sensors are used, a 
parallel development indicates that with increased measurements 


^ u^(z^) for as many of 
n=l 




meas 


i=M+l 


N 


=~I 


n=l 


% 
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yielding an eajiected improvement of 

4 

M r w 


AJ, 




teas 


i^l 


y’"' i 

X ^f(Zn) 

■n=l 




[6.] 


For the second case, of which the mirror protlem is typical, the 
distortion in the hi^er modes is caused hy the control action and the 
measurement error at the n sample point is 


i=N+l 






<ln 


.j=l 



2 

From the properties of the disturbances, the esqoected value of Wg(z^) 
can be determined and is given by 


N r 


E| w^(Zv 


i=l 


"1 2 


nO - 2^ X 


)t 


J1 


J=N+1 




[^64^ 


and 


N U r N r CO 

= I = Is I I 

n=l i=l I n=l j=H+l 


u.(z )i|r.. 

j'' n'^ji 


■[64b], 


This error is on the order of magnitude of the error in the uncontrolled 
modes Jjj and represents the effect of estimating the modal coefficients 
from output measurements. Increasing the number of sensors, as suggested, 
will substantially reduce this error if the remnant Jj^ is dominated by 
the first few higher modes. 



VI. SUMMARY OF lESIGH GOMSIDERATIONS 


The performance index for the system is broken into two parts; 

(l) that contributed by the controlled modes, and (2) that due to the 
uncontrolled modes. Application of active control reduces the error in 
the controlled modes to any arbitrary level while the disturbances 
producing errors in the (uncontrolled/unmonitored) higher modes cannot 
be counteracted. Jhrther, the corrective forces applied by the finite 
number of discrete actuators excites additional errors in the higher 
modes. When mode estimation is employed, errors due to measurement 
imcertainty are introduced into the first N modes. Thus, the total 
system error is given by 


Jl = 


'N, 


+ Jtj + J| 


meas 


E. 



where J^r ^ Jr ; and J-n are defined in Eg,s. I 59a I and 
■^meas ^ "c • ' 

observed that these errors may be minimized' by factors under the 



It was 


designer’s control. These factors are actuator location, pad size, 
and sensor locaition. Selection of actuator location permitted 

. ‘ I ' 

minimization of the excitation of the more significant uncontrolled 

‘ . ■ 

modes. The pad size was seen to act as_a filter which attenuates the 
effect of the control input in exciting the higher modes. The pad 


size is selected to cause the, modal' content of the applied force loading 

* « 

to drop off quickly above the mode. Combined with actuator locations 

^ I » 

that minimally excite the first few (most significant) higher modes, the 
pads together with the plant provide the desired modal filtering for "the 
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remaining high-order modes. Finally, if estimation errors are to he 
minimized, sensor location requirements become the same as those for 
actuator placement. In some instances, additional sensors may be used 
to reduce the errors in estimating the modal coefficients. 

The above points are illustrated most clearly by the example of 
a simply supported, thin, square, flat plate with, the following 
parameters: 

Thickness 0.5 inch (1.27 cm) 

Length 50 inches {j6.2 cm) 

Width ■ 30 inches (76.2 cm) 

Young's modulus 10^ psi (70-5 X 10' gm/cm^) 

Poisson's ratio 0.2 

Table 1 contains values of flc? which are inversely proportional to 
mode transmission as given by the steady values of the A matrix. The 
modes are ranked, from most to least significant, in terms of decreasing 
(l/ti^) and control of the low-order modes is desired. With the modes 
ordered with m on the ordinate and n on the abscissa, contours of equal 
(0:5^)^ become circles in the n,m plane. These circles may be approxi- 
mated by squares^ for exaa^)!©, if 23 modes are to be controlled they 
would include 1 S n — 5 aiid 1 ~ m — 5 as shown in Fig. 17a. The next 
highest mode in each direction is m = 6 and n = 6 which has 25 mutual 
nodes equally spaced over the plate, see Pig. ITti. Actuator placement 
at these nodes (which can always he achieved since the m + 1 mode has 
m nodes) makes orthogonal add precludes excitation of any mode for 
which, m or n equal 6. Alternatively, the first N most significant 
modes may be controlled as’ illustrated' in Fig. 17c « This is achieved 
by the actuator placement of Pig. 17d. The unexcited modes in the 
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317.5 356.7 427.2 536.5 695.4 917.5 1219.7 1621.6 2146.2 2819.4 


536.6 587.2 676.8 


853.9 917.5 1028.7 1194.9 1426.9 1738.5 2146.2 2670.1 3332.9 4160.5 
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Table !•- Values of ixr in 10^ radians'^ /seconds'^ for a flat square plate. 
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Cb) An actuator 
configuration utilizing 
twenty-five actuators. 
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(d) An alternate actuator 
configuration utilizing 
twenty-five actuators. 


me of the controlled^ unexcited 
Lat rectangular plate. 
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latter case are those with m = 8 or n = 8. In any case^ final' 
>de-cennination of the optimal trade-off requires detailed evaluation 
of the performance index. 

Since data regarding the disturbances to which a primary mirror 
STxrface is subjected are not presently available^ a disturbance profile 
characterized by a modal force coefficient with a standard deviation of 

Un , /plate 
_ ^ _ N/area 

‘^^i ~ ~ ~ 586.4 

in pounds per square inch, for all i, was assumed. This profile yields 
an rms figure error of 



F = 


I plate 
area 


w 50 X 10“^ inches 



for the uncontrolled siarface which is in reasonable agreement with the 

figuring errors of the mirror whose diffraction pattern is given in 

Fig. 5* Pad size was selected to be 0.5 inches (1.27 cm) x 0.5 inches 

(1.27 cm). Fig. 18 displays the rms figuring error for a type 1 seivo 

versus the number of actuators for placement of the tsrpe illustrated in 

Fig. ITb. The type 0 error was evaluated for optimal gain but did not 

1 

provide a significant improvement for the error profile considered. 

The details of the procedure used to obtain Fig • I8 'are contained in 


Appendix C. 
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Figure J]_/plate area for a simply supported flat square plate. 
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Tiie preceding sections contain the development and summary of ^ 
design considerations for the discrete control of a distributed- 
parameter system. A simply supported flat rectangular plate has been 
used as an example because it possesses unique properties which clearly 
reveal the results of design decisions which are obscured in most 
distributed system control problems. The following section presents 
the application of the design technique to a plant which is representative 
of a thin deformable mirror and whose con5)lexity is more nearly 
commensurate with that of plants generally encoimtered in practice. 



VII. MODAL CQMTROL OF A FREE CIRCI3LAR PLASE 


Modal Representation 

The equation of motion of the free circular plate of Pig. 19 imder 
forced vibration is 


.-.2 „ d^(r,e_,t) , „ 

V S V w(r,e,t) + p = p(r,e,t) 


[ 68 ] 


where w, p., and are expressed in cylindrical coordinates. Assuming 
solutions separable in r, 6, and t [i.e., w(r,0,t) = Ci(t)fi(r)v 3 _( 0 ')], 
the analysis follows that of the rectangular plate to yield 


where 


d\{t) 2 


dt 


^ + aiLC^(t) = 0 




dV(0) 


d0" 


+ n v.(0.) = 0 


[6910 


r^d^[fi(r)] dfj|,(r) 


dr^ 


+ r — ^ - n^)f^(r) =0 'v [^90 




[ 70 ] 


Since the mode shape given by [6910 is periodic in 0, n is an integer and 
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Ui(r,0) = A(jn(]Sir) + BiIn(kir)]cos (n0 + 0^) [jl] 

where A and 0 are arhitrary constants. J and I are the n"^^ order 
n ^ n n 

Bessel and modified Bessel functions, respectively. The values of 

and Bj^ are determined throu^ substitution of [^Tl] into the boundary 

conditions, which for the free plate, arise from the absence at the 

free edge of both the bending moment in the radial direction and 

( 10 ) 

vertical shearj i.e.. 



V d 
r dr 


n^v 





free edge 



and 


I iJ— + i A 

y,dry^j,2 r dr 


3,2/ ^ r dri r J 


= 0 . 




free edge 


Under these conditions there exists a denumerably inf inite sequence 
of eigenvalues- 


\ = V . W 

for which the associated eigenfunctions form a complete orthogonal set 
permitting, both p(r,0,t) and w(i;,0,t) to be expanded in a uniformly 
convergent series as assumed in Eqs. and §2b]. 

Table 2 contains values of k^^ for several modes. Because of the 
importance of the mode shapes relative to design decisions, the radial 
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components of the first 21 modes, rariked in order of increasing 
frequency, o^, are plotted in Fig. 20. As a result of" the nature of 
the 0 variation in Eq. each pair m,n is associated witli two 

distinct inodes given by 


mode pair m,n = 



(r) cos n9 
(r) sin n9 



for n 0. For n = 0 a single distinct mode exists for each pair of 
m,n. To minimize the expectation of the square surface error, the 
actuators should affect control on the most significant modes as 
determined by the transmission factor X/po^ and the disturbance profile. 
For the purposes of this paper, and as in the case of the rectangular 

p p 

plate, a force distribution with = ag. is assumed such that an 
uncontrolled rms figure error of 50 x 10“^ inches (1.27 x lO”^ cm) 
results. With this white modal disturbance the system objective is to 
exert control -on those modes with the smaller values of kmn Table 2. 
Actuator Size and Placement 

The actuators are again modeled as displacement actuators in 
series with a spring which is attached to the plate by means of a pad. 
The springs are relatively soft to mabs the effects of mirror 
displacement feedback negligible as discussed in Appendix D. The pad 
shape is a portion of an annulus boimded by constant increments in 
radius and angle. The elements hj^^ are evaluated as 











Figure 20.- Modes of free vitratron of a free circular plate. 
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jth 

area 



area 








A[Jn(W> 


+ cos (n9 + ejr dr de 



for which the 0 dependent, component is 


A0. 


sin (nA9:i/2) 
J 

(dA0./2) 


cos + Qq) 


where noimally the increment A0 is constant over all' j permitting "the 

4 

inclusion of this component in the decoupled plant dynamics. The r 
component is somewhat less tractable. Under the substitution , 

n = W 0^0 


the r dependent portion of i3^j "becomes 



+ BIg(Tl)]Tl dTl 


Which is integrabie if n is an even integer^ but requires numerical 
integration or use of tables if n is odd. The effect in either case 
is that the radial component of hj|_j decreases as increases. 
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In determining the placement of the actuators note from Pig. 20 

, f _ • 

that each mode m,n has nodes at m distinct locations along lines, of 
constant 9 and 2n nodes circularly. The placing of 2n actuators 
circularly at equal inteirvals results in their position coinciding 

V 

with the nodes of one of the modes described by Eq,. Consequently^ 

this mode is not excited while the one spatially shifted 90 ^ is. 

To determine the desired actuator placement a knowledge of the 

disturbance profile is required. Under the earlier assumption of a 

white modal disturbance spectrum {cr| = ), the objective is to 

i 0 

control the modes with the smallest values of Inspection of 

Table 2 indicates that lines of constant k ^^ tend to form triangles 
connecting m to n> where 

n ? 2m . [Vt] 

If the controlled modes are n ^ %iax ^ ~ “max controlled area 
of m^n plane is a rectangle (see Pig. 21a) which should approximate 
the region of the first K significant nodes. Actuator placement would 
fall at the nodes of the next highest modes which require 

^actuator ~ ^^°taax 

actuators corresponding to the mutual + l) nodes radially and 

2(,rijjjax l) circularly. The controlled modes are the (2n ^^y + l) 

(“max ■** 'bounded by the rectangle m^^y , plus the m^gy + 1 

excited modes for which n = n^^y + 1. Thus^ the total number of > 
controlled modes is N = 2(n^ + l)‘ (m^ t l) = ^^tuator*' 


The modes 
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not excited "by tliis actuator placement are those at n = n^^^y ^ 

whose nodes fall on the lines of constant 9 where tlie actuators are 

placed (those at n = + 1 spatially shifted 90° are the extra 

. 1 

+ 1 modes included in the N controlled modes). Additionally, 

i 

the mode whose radial nodes are selected as actuator locations is not 

t ) i 

* I ? 

excited. The pertinent controlled and unexcited regions of the m,n * 
plane are illustrated in Fig. 21a along with the tcor re spending' 
actuator placement in Fig. 21h. However, because of the tendency > 
of the lines of constant to form triangles as indicated in Eq,. [ 77 ], 
control of an area in the m,n plane as indicated in Fig. 21c is 
generally desired. This can be accoD5)lished by the actuator placement 
shown in Fig. 21d. Note that in the latter control scheme while the 
N most significant, modes are controlled, it is not possible to preclude 
the excitation of the next most significant modes. The trade-off must 
be made on the basis of actuator location effect on the imaging index 

Jj. 

’While, at the present time, the imaging index has been detailed 
only for the rectangular control scheme illustrated in Fig. 21a, 
preliminary results indicate that control of the N most significant 
modes (e.g., see Pig. 21c) is preferred. The rms figuring error 
F = Jj/plate area based on the rectangular control scheme of Pig. 21a 
is plotted in Fig. 22 versus nvunber of actuators for a plate with the 


following data: 
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(a) The m^n plane for 
the actuator placement of b. 



(c) The m^n plane for the 
actuator placement of d. 



(b) Aq actuator; 
.configuration utilizing 
twenty-four actuators. ‘ 



(d) An alternate actuator 
configuration utilizing 
twenty-four actuator s . 


Figure 21.- The location in the m,n plane of the controlled, unexcited, 
and unmonitored modes of a free circular plate. 
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Figure. 22.- ^J^/piate area for a flat circular free plate. 
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TtLickness , 0.5 inch (1.27 cm) 

Diameter ' 50 inches {76.2 cm) 

Yo\mg's modulus lo"^ psi (70.5 X 10' gm/cm^) 

Poisson * s ratio 0 . 2 

These results indicate satisfactory control, yielding diffraction- 

limited performance, as defined hy an rms figure error of less than 

0.5m. inches, is achieved with less than 2k actuators. This is 

significantly less than the 6 l actuators used in the present laboratory 

(l4), 

model which was determined hy sectioning the mirror into equilateral 
triangles 5.75 inches on a side (the 5*75" x 5-75" X 5.75" x 0.5" thick 
triangles represent a thickness-to-area ratio near that noimally found 
in monolithic telescope mirrors) . 



VIII. COUCLUDUSG BEMARKS 


0?h,e modal expansion technique has "been applied to the prohiem of 

t > 

correcting and maintaining, to the tolerance reguire'd for . diffraction- 
limited performance, the optical figure of a plant representative of 
the primary mirror of an -orbiting astronomical observatory. The modal 
technique has been shown to be particularly appropriate for this problem 
by virtue of its relevance to a useful measure of image quality, its 
ability to decouple the system dynamics permitting simple control 
techniques to be applied, and by the extent of the insight the technique 
■affords into engineering design decisions. 

For distributed plants subject to extremely accurate control, it 
is necessary to consider the effects on system performance of all of 
the modes - not just those which are subject to control. In fact, with 
the error in the modes vinder control reduced to any desired level, the 
major system error was shown to reside in the uncontrolled higher- 
order inodes and this is increased by the control effort applied to the 
lower modes. For this reason the most significant design decisions 
are those related to the effects, of the corrective control forces on 
the higher-order modes. The analysis presented in this paper describes 
the effect of actuator size and location on system perfoannance , factors 
most critical to efficient design. The requisite conditions for 
minimizing the number of discrete control inputs required to achieve 
satisfactory perfoxmance were outlined and then illustrated in two 
design exau^iles. The results for both the rectangular plate and the 
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free circular plate indicate that the thin deformahle mirror can provide 
diffraction-limited performance; further, that this performance can he 
achieved with considerably less actuators than that required' for a 
segmented mirror where the thickness-to-area ratio for each segment 
approaches that normally used in monolithic mirrors. 

The disturbance profile (if data on' the ‘effects of thermal 
gradients, spontaneous release of material stresses, or other factors 
producing distortion of the optical surface become available) can be 
readily Incorporated into the design- procedure. This is achieved by 
using the profile along with the transmission properties of the plant 
(plate) to determine the modal errors and the ST modes, yielding the 
largest errors controlled. Extension to more complex plants (e.g., 
shells), while 3:equiring considerable computing effort, is direct. 
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APEEWDIX A 


Determination of the Eigenfunctions and Eigenvalues of 
a Siniply Supported Plat Rectangular Plate 

The equation of motion of a uniform plate in forced vibration^ is 


o o S^(x,y,t) 

+ p — r-s = 





consider first the homogeneous equation. The modes of free vibration 
will be determined through the separation of variables technique by 
assuming 

w<x,y,t) = w^(t)w2(x,y) {a2] 

Substituting into 


Wj,(t)SV^^2(2C^y) + p ^ 

^t 

dividing both sides by PQWj_(t)w2(x,y) yields 



P W2(x,y) wi(t) ^^2 




Since each side of is a function of different variable s^ both 

p' • - , * 

sides are set equal to cd^ - a positive constant. This yields the 
following equations 


S' 

~P 


^W2(x,y) - cD^W2(x,y) = 0 
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TO 


S w, (t) dwj_(t) 


at^ 


dt‘= 


Rearranging [A5] yields 


a?p Wp(x,y) 

V^wp(x,y) - 0 




or 


(v^ + - ■p‘ij^''2f=^»y) = ° '■ M' 

I 

Tiie solution to [A8] is the sum of the solutions to, each of the products 
of or therefore to 

'[A9]. 

'v^ is the Laplacian in Cartesian coordinates. Equations [A9] and 
are therefore 








V- 




y) = 0 


/ \ 

— o WgCx^y) 

ax^ 



= 0 
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A further assumption is made that 




In this case [All] "becomes 

d^,(x) 


wi^(y) 


dx' 


2 ^3' 


" ~ ^ ' 2 " ' ~ “i\| T = 


0 


Dividing hy w^(x)w]j^(y) and rearranging yields first 


j, ^ , [T . ^ 

^5^^) dx^ dy2 "Ms 





then 


, d^^Cx) I~ . d^j^^(y) 

— + ffl.\ — 7—r 

^5^^) dx^ ^M S W4(y) ^y2 



Since the left side is equal to a function of x, and the right a 
function of y hoth sides must he equal to a constant + (The 

P 2 

choice of sign on is arbitrary since choosing as the constant - p 
will yield the same answers.) Equation [ai5] becomes 


d^5(x) 


dx 


+ 


{■ 



= 0 



d^i^(y) 

dy2 


+ W4(y) = 0 , 



’and 
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The solutions to these equatxons are 


W5(x) = A^sin l\j-[x^± (0 ^^ X +. ij 


[ai8] 


and 


wij.(y) =A2Sin'(ny + Ig),. 


M 


respectively. 

f * * ^ 

From equations ^A12^, Qu.8^, and |~A19^ the solution to equation 


W2(x,y) =AjsinK 


“ ± CDj_ 


M 


X 


+ sin (My + I 2 ) 


where Aj is the product of Aj_ and A 2 * In order to evaluate the constants 
in the boundary conditions must he specified. For a simply 


supported flat plate they are 


(10) 


W 2 (x^y) = 0 for x = O^a 

y = 0 ;h 




% = 0 = - s 






5x2 


+ V 


ay2 


[a^ 


x=0 

x=a 


and 
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0 = - s 


' 3^V2(x,y) 




ay2 


+ V 


^^W2(x,y) 


y=0' 

y=b 



in light of [a 2]^ the equations Qa 20 and [a 2^ may he written as 


^%(x,y) 


3x2 


= 0 


1:^24] 


x=0 

x=a 


and 


^ WgCx^y) 




= 0 


[a25] 


y=o 

y=b 


Using the conditions for x = 0, y = 0 from equation [a 21^ in [a2o 3 
yields 


ll ^2 ^ 


(a2^ 


Using the condition of equation [A21] for y = h yields 


sin |J.h = 0 


[a2t] 


and therefore 


[a — -:g- n - 0^ 1^ 2^ ... 


[a28] 



From the condition for x = a 


sin 



= 0 



or 



“ = [a30] 


Substituting [A28^ and ||a 36J into [^<0 gives in [a3i] the shape of the 

f 

modes of free vibration^ or eigeirfunctions^ of the plate 

=A^sin5^ sin n,m ~ 1, 2, 3, . .. , [a5:Q 


The resulting mode shape is identically zero for either m or n zerO; 
consequently, equation |A53^ is valid for the range of m and n 
indicated. 

That equation satisfies and [A25^ may be verified through 

substitution. Further, may be chosen to satisfy the requirement 


r f 0D?(x,y)(3x dy = 1 
Jq Jq 1 


[a52] 



ab 



which yields 
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^ 2 (^.y) = ^ sin sin 


rn n S '^^2(2c,y) 

Since from [a^ o| = - 


4 = 


a2 S2 




3x^ 


V 


w 


2(2:^y) 


w, 


2(^^y) 





This could also have "been obtained from [a 27 ], [a 28 ], and |a3o]. 
Since the solution to (6) is 


Wj^(t) = COSTCO A + . [a3^ 

The most general solution to the' homogeneous form of [ai] is 

t. t ^ 

00 

/ N V / X . ^ -t 

w(x^y^t) = > A cos(ob + ^ )sin sin — — \P39\ 

\ JO? j mi ^ mn ^imi' a , Id 

iri;n=l . ’ ^ ^ , 

I * ‘ ' • t 

Thus the eigenfrmctions or modes of free vibration and their 
associated eigenvalues have been determined for a simply • supported, 
flat, rectangular plate and are given by equations ,[ a 34-^ and (A37^ 
respectively. 



apeehjdik: b 


Determination of a Set of Actuator Locations for Which For 
a Simply Supported Flat Rectangular Plate is .Oirthogonal 

The purpose of this appendi:x is to indicate^ for a simply supported; 

flat; rectangular plate a set of actuator locations whidh make the 

matrix (eq. H) orthogonal. The equation of motion ' of , a heam^ is 


Jr p 

2 ^ w(x,t) ^ ^ yw(x,t) 
^ dt^ 



The method of separation of variables is used, consequently; 

w(x;t) = w^(x)w 2 (t) Qb2] 

Substituting r?=] into |bi^ gives, after dividing both sides by 'w(x;t); 


dS/gCx) ^ d^WgCt) 


■"iC^!:) dx^ W2(t) 


dt'^ 


[B3] 


Both sides are set equal to a constant 7^ 


d^Wi (x) , 

- 7^wj^(x) = 0 

dx^ 

d%{t) , 

— + 7 ■W2('fc) = 0 

dt^ 


[^b4^ 


M 


The solution to M xs 




j(t) = sin T^t + Cg cos 7 ^t 


B5] 
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The solution to the remaining equation, will yield the modes of 
free vibration of the beam. Since the exact solution will be used in 
the sequel this information will be obtained first. 

Equation (b 40 is factored as 






and the solution is the sum of the solutions to 


(5 ■ ' ° 



and 



+ — Pl(^) = ° 



Thus 


w-j_(x) = Cjcos X + Ci|. sin +0^ cosh +C 5 sinh 

at this point the boundary conditions are brought in. Eor a simply 

( 10 ) 

supported beam 


w-j^Co) = 0 

(bio^ 

11 

0 

|bio0 
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and there are no moments at either x = 0 or x = b 


dx^ 


= 0 


x =0 

x=b 


These conditions require 


0 = 05+05 




0 = C^cos — ^ h + Ck sin b + Cc cosh b + CfCsinh b [bII^ 

^\pF ^ ° 

[bIIo] 


T 


0 = - ^ 03 + ^ C 5 


0 = - 


if 

fF\ 


C2 cos — ^ b + Ck sin — ^ b 

^ ff ^ \f¥ 


y2 

'W 


r 


Cccosh — ^ b + sihh — ^ b 

^ \[T ^ \|T 


From equations [b 11^ and [b 110 it is determined that 




= C5 = 0 


[bis] 


leaving 


0 


= Ck sin — ^ + C< sinh 

^ \TT ^ 


-2Lb 

\T? 



7b 

0 = - C ij. sin + Cg sirih 


\P 


\F 



A nontrivial result requires the determinant of the coefficients of 
O 4 and Cg in equations [bI^^^ and [b 1^ to be zero. Consequently; 



Substitution of [bi 6] into [b 15]] indicates that Cg = 0 and 


7 _ nut 
\/T " b 

The eigenfunction of the homogeneous simply supported beam is 



w-j_(x) = Cl). sin •— X 



The solution to the problem is now considered through the method 
of finite differences. A number of stations are located at eq.ua! 
intervals, Z, along the beam as shown in the following sketch 

2 ~ 5 | 
t J. 

0 xx X2 Xji-i Xn 

and the equation 



d^wW 





is written at each point, where w® = col cd(xx) ^ • ®(Xq). The fourth 

derivative approximation used can be obtained by first obtaining the 
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Taylor series expansions ai)out x^ of 'w(x^^2.)^ ^(^+2)^ ^ 

w(xj._2)* These exj)ressions are 

w(x^l) = w(x ) + w’(x )Z + w"(x )— + . . , + w^^\x )— [b2o5 
r^x r r r ^ n! 


w(Xr-i) = w(x ) - v'Cx )Z + w"(:x )i~ - . , . + 

2! , , n! 

. ' [b20^ 

♦ t 

w(x 3;.+2) = ■w(^) + tr’(xr)2Z + w"(xj.)“^ . . . + 

i > 

' • ■ ' [b2O0 

r ^ ^ 

.2 ^ ‘ ' n 

. . ' (b2o5 

subtracting four times (|^20a^ + j^20b^)* fromp20(^ + ,|B20<^ yields 

^ (2^ . It) 


- ‘^ »(Xi.+l) + * "(='r+2> * ■ ^ 




+ . » . t 2 


d^(x) 


ax' 


,2m 


(2^“ X) 


(2m)! 




+ . . . 


[b2i] 


After rearranging equation [b2i] "becomes 



8l 




dsT 


+ « . . 


X=3C^ 


+ 2 


d^r(x) 


dx‘ 


2m 


^2(m-2)22(m-l) 

i^l 


+ . . . 


x=xr 


[b22] 


The term on the left will he denoted The finite difference method 

approximates the plant equation as 




T 


= '0 



Collecting the expressions for each point and arranging them as 
indicated hy the definition of w®^ yields the matrix equation. 



%dm • 





w(xi) 

w(Xn) 



The solutions to this equation are the eigenvectors of the finite 
difference - matrix representation of the plant, matrix 

*■ ‘ * i * 

is a real symmetric matrix and consequently, has orthogonal. 


eigenfunctions. 

' ' - r • '' ' 

If the right-hand side 'of equation jB22j is used instead of-tihe 

, * * * 

left it is possible to obtain' an kialytic expression for the numerically 
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obtained finite difference answers. That is, consider the equation 




0 = ^ _ 


cbc^ 


w(x) + . . . + 2 
T ^2di (2m) ! 


d^w(x) ^(m-1) ^ 



The answer to this equation evaluated at x = Is equal to the finite 
difference answer. 

Consider the function u ^ = Ci^ sin x which, from [bi 8^, is seen 
to he the exact eigenfunction. Substituting this expression into []b 2^ 
yields, after cancelling u^ 


0 = 


2^ + 

t2 ” t2 


+ (- 1 ) 



z2(m-2) 

(2m)! 


22(m-l) 


• • • 



Thus T3i(x) is the solution to [b 2^,, and u^(xx) = 'w(xj_), when 




(2i)2(m-l) ^ 


m ^ 3 



I is the separation between adjacent stations on the beam and it can 
be made as small as desired. Consequently, equation [b 2'^ indicates 
the resultant convergence of - — . to More significant for the 

purpose of this appendix is the fact that the i component of the 
eigenvector of the finite difference solution is equivalent to 

^ t f 

the eigenfunction of the beam evaluated at^a location' corresponding 
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to the point in the finite difference representation of the heam. 
Since the matrix of Eg. is 


U]_(X2_) . . . Ujj(x3_) 





each column corresponds to one of the orthogonal eigenvectors of 
and the matrix is itself orthogonal. 

To this point the proof has conceraed the solution for a beam, 
while it is desired to show ^ 

^ I, 

N i .... . ' 

^ U^(X35.)Uj(Xjj.) = 0 ‘ ^ 

E=1 . ‘‘ ' 

for the plate. For the plate, Appendix' A shows that' 


/ ^ * I . mi 11%: . * 

u.; (x^) = Ca Sin — — r- sin — - — 



conseguently. Eg. may be^ written 


N K 

V Ui(xj^)uj(xj^) =V 
k=l k=l 
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Consider an array of locations 


JSf = XY 



with. X locations in the x direction and Y in the y direction. Thus 
hecomes 



X , ^ , 

A^=!L Bi^=l 


Each component in the second expression for the right-hand side of (b35^ 

, t 

xs equivalent to the heam^ consequently^ either the first or second 

t 

summation will he zero unless i - J. 

Equally spaced points will provide an orthogonal matrix for the 
purpose of relating performance specification in the original and 
transformed systems. 



APPENDIX C 


Evaluation, of the BMS Figuring Error for a Simply 
Supported Flat Rectangular Plate 

The purpose of this appendix is to describe in detail one of the 
procedures used in determining the results contained in Figure l8. As 
indicated in equations Cte] and [67] the exact determination of Jj would 
require the evaluation of an infinite number of terms. In using a 
finite number of terms p to approximate Jj and, therefore the rms error, 
it is desired to select such a P, if possible, which would place a bound 
on the amount by which the approximate value of "Jj would differ from 
the true value. 

The value of P which should be used is a function of the individual 

t 

terms in the sequences 

ai, ag, • • (C]] 


and 


Since 


2 2 





one procedure might be, for monotonically decreasing' values' of ‘c^gg, 

to observe values of Oigg in for increasing i until values of 

C-, are obtained which are significantly less than the accuracy desired. 
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The value of Jj may be determined for two values of P in this neighbor- 
hood to determine whether or not has been obtained accurately enough. 

While such a procedure would be adequate the specific nature of 
the present problem permits the selection of P on a more rigorous 
basis. For a square plate the expression contained in Appendix A for 
the eigenvalues becomes 

[eg 


r „2 . „2 

j — m. +11. 

^ pa^ L^' iJ 


Prom Eqs. ^2^^ 


r So ^ 

Since m and n each take on the values of all the positive integers 
the right-hand side of [c 53' may be rewritten , 

f * ‘ » 

* J 

CO CO CO 2 

In order to remove one of the infinite sirmmations , use is made of the 
symmetry of the eigenvalues by writing [C6] as 


Vi 
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where is considered a constant as has already "been assrmed (see 
Eq. [6^). The terms i^ + in the series 


II 

■ 1=1 


Q.2 + 




are placed into correspondence with the positive integers in the order 
indicated by [c^. It is desired to indicate for the i and j 
corresponding to the value in [ 08 ] that 


In the first n values of i there are 


■n 


i=l 


terms where T is 


rp“= 

• ' 2 




? 2 

The j = 1 term is the miniinum value of 1 + j for any value of i. 

* * , , 

Since, for i = n ’ 


n^ + 0^ > n^ 



if 


n2 



2 
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inequaXity 09^ will Xe obtained. Since is valid for 


n > 1 


then [c^ holds and the right side of [C7] may be. written 


[eiq 


[sJ 


'^ < y y ^ < 2 !-!iiy x roiq 


The advantage of the last series is that its sum has a bound. The 
partial stmi for the first 2^-1 terms is 


^■t-1 - 


1 . 




+ 


1 

(2^ - 1)^ 


which is less term by term than 








1 


( 





which is 


thus 



Q' 


2‘fc-l 


1 



1 - 


1 

25 


25 

25 . 1 



25 - 1 
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The second of the two terms is always negative and goes to zero as 

* 

2 '^ - 1 increases. Thus the siim of the primed series from the 2 “ - 1 
term to the end, Q , is less than the second term on the right side 

of 0-9] 


«K< 



2^-1 



Since the primed series is greater than by term than the k series, the 
sum of the k series over the same terms must also he bound by the same 
amount. 


CO 


I 

k=2-fc 




The original series is smaller than the k series teim by term, 
therefore, since 



w2(x,y,t)dx dy 


2 -1 


=E 


^ 1 


i g + nf] ^ 


i=:l 




then 



2 -1 

w^(x,y,t)dx dy < ^ 
i=l 




25 - 1 
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and 


n "b 

(x,y,t)dx 


2^-1 


a^a? 

1 


8 


a? i 


t-1 


< a 


1=1 K ■*• ^Pt ^ 


[c2ir) 


The inequality may be used either to place a bound on the error' 
for a fixed number of modes considered in the evaluation^ or conversely, 
it may be used to select the number of modes which must be considered 
to keep the error below a certain level. It should be noted that the 
value of in [oaq is the result after control and must be chosen 
conservatively enough to reflect the induced error as well as that 
already present. Once the bound on the error has been selected t may 
be determined from |^2i0. This equation was based on the use of 2^ - 1 
terms in the series, consequently a square n x n array of modes cannot 
be used for which 

SiEiil > 2* - 1 
2 

Finally} the tightness of the bound is dependent on how closely ' 
the k series approximates the true series for the terms after the 2^-1 
term. The bound may be tightened by observing actual deviations and 
adjusting the inequality [c 24] by the appropriate amount. 



APPEHDIX D 


Mirror Displacement Feedback 

Tlie actuators considered in this appendix are modeled as types 
which are being considered for actual usage-^^"^^. This model consists 
of a pure displacement actuator acting aga'inst a spring and a backing 
plate which is stiff relative to the spring. To obtain a specified 
force the displacement actuator is commanded to' a new location relative 
to the undeformed mirror. 

If the mirror has deformed, the displacement of the mirror will alter 
the magnitude of the applied force. It is assumed, arbitrarily, that 
one form of this displacement feedback might take is 

Pi(x^y:>t) = a^(t)p^(x,y) -Kw^(x,y,t) j^Dl] 

where Pj^(x,y,t) is the force density applied by the i^^ actuator, and 
where 

Wi(x,y,t) = w(x,y,t) (^DsJ 

oVer' the i^^ pad 'area , and 

* * I 

‘ WiCx^y.t) = 0 [d210 

V 

elsewhere , and 


a^(t) 


KZ.(t) 


If 


y)dx dy 
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w^(x,y,t) may te e 2 q)anded as 


w^(x,y,t) 


=I 

j-i 


cij(t)u3(x,y) 



where 



w^(x,y,t)u^(x,y)dx dy 



Since (x^y^t) is> zero except over the pad area of the i'^li pad I D5 j 
can be written 


oij(t) 


-If w(x,yjt)n^(x,y)dx dy 


pad 

area 



Substituting the modal expansion of w(xjy,t) yields 

* CO 

J f Yj M 

1'^'^ pad ^ ^ 
area 


Interchanging the order of summation and integration yields 


00 

'ij f I 


k:=I 


itb pad 
area 
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5}he total feedback force, fg(x,y,.t) is 


fg(x,y,t) 




Substituting [d^] into []d0 


i=l 0=1 


Interchanging the order of summation 


CO n 

^ X Z 

0=1 x=l 



By comparison with equation [l2b] 


nan be expressed modally as 


fs(x,y,t) = X (t)iij(x,y) [i)12] 

d=l 


where 


N 

a!(t), = - K X clj(t) [d13] 

i=l 


Equations [d 133 used to put the expression for aj(t) 

into a different form 




N c» 


"XX f / 

i=lk=l y- 


i'^k pa.^ 

area 



K 



9k 


IntercTianging the order of summation 


K 


=^Ck(t)^ J J Uk(x,y)uj(x,y)dx dy [pi^ 

k=l i=l ith pad 
aref 


thus 


- a'{t) = +KZ c(t) 

where Z is an to X CO matrix which has individual elements of 

N 

i-l ith pad 

area 

and where the spring constants are all assumed equal* 

The matrix Z is located in a local feedback loop around the 
diagonal plant matrix since it describes amplitudes in the force modes 
as a function of the displacement mode aii5)litudes. Since^ in general, 
the elements of Z are non-zero, tMs feedback causes the system to 
become coupled. In this particular application the coupling exists hut 
can he made to have negligible amplitude. Since the displacement of 

* t 

the mirror is on the order of microinches the spring constant is ad-justed 

• I * i 

so that the actuator throw reg.uired is on the order of Inches. 

Alternatively, a form of spring feedback may occur whidh can be 
treated without introducing coupling. Suppose that the expression for 
Pi(x,y,t) is 
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Pi(x,y,t) 


K(idi “ imjpifcy) 



3^(x,y)djc dy 



which would he obtained by decreasing the pad area until the mirror 
displacement over the pad area is constant. This, depending on the 
manner in which the pad is bonded to the mirror, appears to be a 
reasonable assumption . Expanding [di 8] yields 


Pi(x;>y^^) = 


3i(x,y^ 


Klra^3i(x,y) 


J/Y^i(x.y)te dy fj 3^(x,y)dx dy 


[d19]' 


This e( 3 uation can be analyzed by a procedure similar to the preceding 
paragraph, or equivalent results may be determined from an inspection 
of the appropriate block diagram. Proceeding as previously, the second 
term in [^ 19 ^ is 


i“sj,(x^y^t> 



= - Klu^%(x,y)/ J J ^i(x,y)dx dy 



where is the constant value of w(x,y,t) over the pad area. It is 
desired to express fg(x,y,t) in a mod^ expansion 


where 


= ^ayt)uj(x,y) 
. ‘5=1 




w.(x,y,t)p (x,y)u ('x,y)dx-dy 

[D22] 



96 


Since Wj^(x,y;t) is constant over the pad area and' is zero 

elsewhere [d 22] ■becomes 




- Klin^(t) 


ff ^i(x,y)'lx ^^Srel^ 


f J ^i( 2 ^^Huj(x^y)dx dy £t>2^ 


The integral has been previously evaluated as h j (El* ) } therefore 


, , - KZffi.(t) 

aj .(t) = h 


J J" ^^(x,y)dx dy 


•Ji 


[D2g 


and if the integral of p.* (x^y) is equal to that of Pi(x,y) then 

X J 


a^(t) = 


- K 




m\t) 


y)dx dy 



where 


2^t) 


w(xj,y^,t) 


= w^^(t) 



y(xj^yw^t)J' 


In this case the general result indicates that there is coupling in the 
system. If the system is ass\imed to have only finite eigenfunction 


content then 
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N 


y) 


[d27] 


i=l 


If [d 27] is placed in matrix fonn 

c(t) 


[ d 28] 


vhere 


# = 


Ui(xi,yi) . . . %(xi,yi) 




[d29] 


The spring feedback loop as described by j^l^ and |D2^ is shown in 
Figure D-1 



Figure D-1 * 

The spring feedback loop. 





98 


As indicated previously^ under appropriable, conditions^ the matrix can 

^ * * \ 

be written as the product of a diagonal matrix and a riondiagonal matrix 

\ ‘ * i* 

as in 34 • In this case the spring feedback loop becomes that shown in 

- t 

Figure D-2 



Figure D-2 


Tb.e spring feedback loop for h^=gV. 


Both, 6 and 


are diagonal elements. Consequently, 


P^(x,y)dx dy 


when ^ the system is. decoupled (see Fig. D-3) and the effects 


of the spring feedback can he readily included In the system design. 



Figure D-5 

The decoupled spring feedback loop. 
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In stumnary, under a specific set of assumptions, the effects of 
the mirror displacement feedback can he treated without introducing 
coupling effects into the system. In general, coupling effects are 
present, however, the mirror displacement feedback is rendered negligible 
through an appropriate choice of spring constant. 
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